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1. PROBLEMS AND PROGRAM 

The finite temperature transition of QCD can be 
seen as a change in the structure of the hadrons 
and as a symmetry breaking transition - a change 
in the structure of the vacuum (we shall take the 
most economical attitude that deconfming and 
chiral symmetry restoration are related). These 
phenomena are observed differently and carry 
complementary information. We aim at a corre- 
lated analysis involving hadronic correlators and 
the vacuum structure including field and density 
correlations, both non-trivial questions. 

To understand the hadronic phenomenology at 
T > (see, e.g.0) we need to describe the dom- 
inant low energy structure in each channel. We 
must be prepared to cope with: wide structure 
replacing the well defined, T = pole; the dif- 
ficulty of separating this structure from the rest 
of the spectrum; non-isotropic dispersion law - 
all intrinsic (and relevant) physical aspects. Due 
to breaking of the Lorentz invariance at T > 
we must study the general correlators. However, 
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in lattice calculations the finite time extension 
l T = 1 /T prevents the disappearance of the higher 
excitations in the time propagation and more re- 
fined analyses are required. Likewise, for the de- 
scription of the vacuum structure one must disen- 
tangle the physically relevant structure from UV 
fluctuations. We use a "gold - washing" cooling 
algorithm with nearly scale invariant instantons 
above a short range cut-off po — 2.3a M. Since 
instanton - anti-instanton (IA) pairs annihilate 
in any cooling their study is more sophisticated. 

To allow for high T with large N T but moderate 
N a and (3 we use anisotropic lattices: 
a a /a T = i>l, T={i/N a )a- 1 (1) 
This ensures a fine discretization of the time axis 
and thus more detailed information. It also allows 
a fine variation of the temperature at fixed /3. 

We approach these questions first in quenched 
QCD, which shows a deconfming transition. But 
since chirality directly involves the quarks, in 
their absence the temperature effects on the 
hadronic correlators and on the topological struc- 
ture may be delayed or modified. Therefore the 
final aim must be a full QCD analysis. 
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2. ANISOTROPIC LATTICE STUDIES 

Calibration and scaling: Anisotropic lattices 
are introduced with help of a (bare) coupling 
anisotropy 7, i.e. for QCD with plaquette action: 



Sym = — 



i?eTr(P (T(T ) + 7i?eTr(P CT 



(2) 



with P M „ = WV(ll) 



Euclidean symmetry 
should be recovered for physical quantities when 
expressed using the cut off anisotropy £ ||,|| : 
FZ(z) = FZ(t = £z) (3) 
This fixes £ given 7 ( "calibration" ) . If the action 
leads to strong artifacts, this relation cannot be 
fulfilled simultaneously for all observables. In a 
first analysis of these effects we use the tree level 
improved actions defined with P in Eq.(2) given 
now by the sum of loops: 

Pp, = coWV(ll) + cxW^U) + c 2 WV(22) (4) 
with the 14 7 (12) loop averaged over directions (it 
is easy to see that at tree level the bare anisotropy 
affects all loops similarly). We compare (1): Wil- 
son action, (2): Liischer - Weisz Symanzik action 
c = 5/3, ci = —1/12 and (3): the "square" 
Symanzik action cq — 16/9, c\ = —1/9, c 2 = 
1/144 f§. In Table I we present SU(2) results on 
8 3 x 24 lattices at 7 = 3 and (3 = 2.339, 1.768 
and 1.772 (4000, 4000 and 2000 configurations 
respectively, separated by 10 sweeps after 10000 
thermalization sweeps; x 2 cannot be compared 
between the different actions). The parameters 
are chosen such as to have the same cut off do- 
corresponding to (3 = 2.25 for the Wilson action 
at 7 = 1. (Notice that for the Wilson action 
A(7 = 3) ~ 0.8A( 7 = 1) §.) We fit Eq. (3) 
choosing for F^{m^) planar Wilson loop ratios 

R nA m n) = W <?p,( n <7: m p,)/ W >?li( n o- - (5) 



for 



= 1, . . . , + 1, fj, = a,T. For the 



Wilson action £, 



pert. 



3.3, hence we have rather 



large non-perturbative corrections. The tree level 
improved actions already seem to reduce both 
the non-perturbative effects and the lattice ar- 
tifacts. Results for SU(3) and for non-planar 
loops and physical isotropy checks for instantons 
on anisotropic lattices will be reported elsewhere. 
Deconfining transition: Once the calibration has 
been performed at T = it is assumed that £ 
does not depend on the lattice size, and therefore 
we can increase the temperature by reducing N T . 



action 


ifc 




i?4 


i?. 


(1): € 


4.02(6) 


3.90(2) 


3.89(2) 


3.93(4) 


X 2 /d..f. 


86 


2 


.4 


.2 


(2): £ 


3.52(3) 


3.44(1) 


3.43(2) 


3.42(2) 


X 2 /d.f. 


33 


1 


.3 


.2 


(3): £ 


3.51(1) 


3.46(2) 


3.46(3) 


3.45(4) 


X 2 /d..f. 


36 


10 


3 


1.3 



Table I: Cut off anisotropy for 5/7(2). 

In Fig. 1 we show the Polyakov loop suscep- 
tibility as function of 7 for 3 lattice lengthes 
N T = 20, 18, 16 at f3 = 5.68 (pure SU(3) theory 
with Wilson action, N a = 8). 




Figure 1. Polyakov loop susceptibility \ vs 7- 



Mesonic correlators near T c : Methods for pro- 
jecting onto the low energy part must be carefully 
defined, since a simple reweighting of the spectral 
function may deform a wide peak (while it would 
not change a pole). Using the quark propagators 
S we measure the "wave function" 
F r (P,:r,i)=y>' Pz "(yi.ya)* 
yi,y2 



(Tr [TS(yi, 0; z, t)TS(y 2 , 0; z + x, t)}) (6) 
Here T defines the channel (it, p) and w(yi,y2) 
is the source, to be determined iteratively, aimed 
at optimizing the signal of the ground state. Pure 
SU{3), 12 3 x N T lattices at /3 = 5.68 and 7 = 4 
are used (Wilson action). At this 7 one finds 
T ~ 0.93T C , 1.03T C and 1.15T C with N T = 20, 18 
and 16, respectively (see Fig. 1). The calibra- 
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tion has been done with N T = 72 (T ~ 0) and 
found £ ~ 5.9 from Wilson loops. On this T ~ 
lattice we measured also quenched pion propaga- 
tors in space and time directions. They are found 
to show the same cut off anisotropy £ = 5.9 for 
7i? = k t I ' K a = 5.4 in the fermionic action for Wil- 
son quarks (k: the hopping parameter). About 
20 configurations at N T = 20 and 18 have been 
analyzed, N T = 16 is under way. We pursue a 
number of strategies: 

(a) Using F r (0,x,t) from a simple source like 
point ("pp": w(yi,y 2 ) = S yi0 S V2 o) or wall ("W: 
w(yi,y2) = 1) we fit an ansatz with three poles 
corresponding to the ground state exp(— ax p ) and 
two radial excitations || . We obtain in this way 
the wave function parameters a, p and a first es- 
timation of the lowest mass. These a, p are then 
used to project onto the ground state at the sink. 

(b) Using the same a, pa new, shell model type 
source is constructed by smearing one ("ep") or 
both ("ee") quark propagators with exp(— ax p ). 

(c) "Effective" masses are extracted by fitting a 
cosh around each t for F r (0, 0, t) for the various 
sources and sinks. 

(d) We make an analysis of F r (0,0, t) corre- 
sponding to binning of the spectral function, 
again using the various sources and sinks and 
checking the stability of the low energy structure. 

(e) We analyze -F r (P, 0, t) for the dispersion law. 
This program is presently in work. Partial re- 
sults from analyses on smaller lattices have been 
presented before ||. Now we show in Fig. 2 
the effective mass plots for various sources and 
the wave functions F r (0, x, t) for the "ep" source. 
There is a strong dependence of the effective mass 
on the sources, even in the region where it seems 
to saturate. This signals a low energy structure 
of significant width. A second observation is that 
there seems to be little change, both in the ef- 
fective mass and in the wave function, inside few 
percents around T c . Typical wave function pa- 
rameters are a = 0.4, p = 1.35. However, be- 
fore interpreting these results we want to do the 
full analysis, also of further quantities such as the 
scalar propagator and the condensate, and at a 
higher temperature N T = 16. Also a signal from 
free quark propagation, again a wide structure in 
the spectrum, may show up above T c . 
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Figure 2. Pseudo-scalar effective mass vs t (upper 
plots: "pp" circles, "ep" diamonds, "ee" squares, 
"ww" crosses) and "ep" wave function vs x, t = 
1,3,5,7,10(9) (diamonds, bars, squares, crosses, 
hexagons: lower plots), at N T — 20 (0.93T C ) and 
N T = 18 (1.03T C ). No reliable estimation for 
the effective mass errors is available; they may 
be compatible with the "ep" - "ee" difference. 
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